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Abstract 

SAR  Interferometry  explores  the  height  structure  inside  a  Radar  pixel.  Two  antennas  separated  with  a  so 
called  baseline  distance  B  in  cross  velocity  direction  take  two  independent  images  from  the  same  scene  and 
compare  the  two  phases  which  result  from  the  same  pixel  at  the  two  antennas.  This  phase  differences  com¬ 
pare  to  different  distances  between  the  pixel  and  the  two  antennas  and  from  this  difference  the  height  of  the 
respective  pixel  in  relation  to  the  altitude  of  the  radar  over  a  defined  level  can  be  determined. 

The  Interferogramme  formation  will  be  explained  and  the  respective  terminology  will  be  presented.  In  this 
context  basic  relations  and  important  geometrical  and  electrodynamics  system  parameters  will  be  presented 
and  defined,  the  phase  unwrapping  will  be  considered  as  well  as  basic  interferometry  processing  procedures 
and  image  formation  principles.  Coherence  is  the  basis  for  Interferometry  as  well  as  a  careful  calibration. 
Therefore,  the  role  of  interferometric  coherence  for  image  interpretation  will  be  shown.  Errors  caused  by 
phase  noise,  de-correlation,  co-registration  misalignment,  atmospheric  disturbances,  baseline  limits,  etc  will 
be  considered.  They,  principally,  are  the  main  reasons  for  coherence  degradation.  Different  Interferometry 
modes  Cross  Track  and  Along  Track  Interferometry,  single  and  dual  pass  Interferometry,  Differential  Inter¬ 
ferometry,  and  Permanent  Scattering  Interferometry  SAR  will  be  discussed. 

1.  Introduction 

Goal  of  this  lecture  is  to  present  the  basics  of  SAR  Interferometry  in  order  to  prepare  the  bottom  for  the  fol¬ 
lowing  lectures  on  Polarimetric  Interferometry  and  its  applications.  An  Electromagnetic  Field  is  character¬ 
ized  by  Amplitude  E0,  Frequency  f,  Phase  cp  and  Polarization  p,  the  latter  describes  the  Vector  character  of 
the  field  E=(Ex;Ey)  principally  as  functions  of  time  t.  Ex  and  EY  are  orthogonal  to  each  other  and  character¬ 
ize  the  polarization  base  of  a  Transversal  Electromagnetic  wave  (TEM  wave)  which  has  no  component  in 
propagation  direction.  It  is  sufficient  for  a  principal  understanding  to  consider  in  a  first  step  “scalar”  Elec¬ 
tromagnetic  Waves  and  there  behavior.  Tat  means  without  limitation  of  generality  one  component  of  the 
Electric  Field  Vector  only  will  be  used,  and  polarization  effects,  principally,  will  be  excluded.  The  scalar 
Electromagnetic  Wave  is  written  as: 


E  =  E0ex(i(27if  +cp))  (1) 

2.  SAR  Interferometry  Principles 

SAR  Interferometry  is  an  established  technique  for  rapid  and  accurate  collection  of  topographic  data,  which 
is  essential  for  establishing  Digital  Elevation  Models  (DEM’s),  it  allows  the  vertical  location  of  the  effective 
scattering  centre  inside  the  two  dimensional  horizontally  oriented  resolution  cell. 

Normally,  SAR  Interferometry  works  with  two  antennas  separated  in  cross  velocity  direction  characterized 
with  the  baseline  distance  B  taking  two  independent  images  from  the  scene  observed  and  comparing  the  two 
phases  which  result  from  the  same  pixel.  This  offers  the  capability  to  generate  topographic  maps.  Adding  a 
third  flight  at  the  same  passes  as  before  and  subtracting  the  two  results  from  each  other  is  called  Differential 
Interferometry. 

Due  to  the  distance  between  the  antennas  the  distance  to  a  single  point  scatterer  on  the  ground  is  slightly 
different  by  the  amount  of  r.  This  slant  range  difference  r  corresponds  to  a  phase  difference  ®  comparing  the 
received  signals  of  both  antennas.  After  detection  of  ®  the  local  terrain  elevation  above  a  well  known  refer¬ 
ence  area,  a  geoid  for  instance,  can  be  determined  according  to  the  following  equations,  which  are  based  on 
the  geometry  for  SAR  Interferometry  depicted  in  Fig.  1. 
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Fig.  1  SAR  Interferometry  Principle 

In  practice  two  different  methods  are  applied  to  obtain  an  interferometric  baseline: 

•  One  pass  Interferometry:  In  this  case  both  antennas  are  mounted  rigidly  to  the  same  platform  with  a  fixed 

baseline.  There  is  one  transmitter  for  the  two  receiving  channels.  An  interferometric  image  is  obtained  af¬ 
ter  one  pass  of  the  platform.  For  the  Phase  difference  results  from  the  two  phases  cpi  and  cp2: 

a  ~  2i?,  +r 

(Px=AKX  ’  <Pl  =l7T  A 

®  ipass  =  (p2-(Px  =  2^Y 

•  Two  pass  Interferometry:  Only  one  Transmit-Receive  channel  with  one  antenna  is  used,  but  the  platform 

has  to  pass  the  same  terrain  twice  in  the  same  direction  along  a  trajectory  parallel  to  the  first  one.  In  this 
case  results  for  the  Phases: 
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For  the  Phase  difference  Onj  will  be  written  in  the  following: 

r 

0  =  m  2  n  -ij- 

m  A 


(2) 


With  m=l  for  the  one  pass  case  and  m=2  for  the  two  pass  case.  The  propagation  path  difference  2 r  is  usually 
much  grater  than  the  wave  length,  therefore, O  normally  shows  an  ambiguity  of  many  cycles.  In  any  case  it 
can  be  determined  mod  2n  only,  the  phase  is  wrapped. 
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For  topography  purposes  an  across-track  arrangement  of  the  two  sensors  is  essential.  Flowever,  along  track 
arrangements,  principally,  allow  the  measurements  of  the  ground  velocity  if  targets  ore  pixels  respectively 
move  between  the  two  different  imaging  times  between  the  two  channels.  Both  modes  are  depicted  in  Fig.  2. 
The  coherence  between  the  two  signals  is  essential  in  any  case. 


Fig.  2  Geometry  of  Across  -  and  Along  -Track  SAR  Interferometry  Modes 


2.1  Basic  Interferogramme  Relations 

From  the  geometry  depicted  in  Fig.  1  the  following  basic  Interferometry  relations  result: 

z=H -Rcos  $ 


After  some  manipulations  results  for  the  slant  range: 


The  imaging  Geometry  is  given  by 


R= 


2Bsin{3-%)-2r 


sin  ( 


s-4)= 


2BR 


r 

B 


(3) 


(4) 

(5) 


The  ground  range  is  determined  by 


y =R  sinR 


(6) 


The  local  terrain  height  z  above  the  reference  area  is  therefore: 


z=H- 


2B  sin(3-%)-2r 


cos  3 


(V) 


In  connection  with  (2)  equation  (7)  delivers  the  desired  Interferogramme  as  function  of  the  look  angle  9.  In 
connection  with  (6)  equation  (7)  delivers  the  ground  range  representation.  The  propagation  path  difference  r 
for  the  both  signals  at  the  two  antennas  is: 
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r=m 


R2  +B2 -2BRsin(  3-£)-R 


L  J  (8) 

m=lfor  one  pass  Interferometry,  m=2  for  two  pass  interferometry  with  separate  parallel  flight  passes.  The 
corresponding  phase  difference  O  is: 


_  .  r  m  2  n 

0=m2 n — = - 

A  /l 


R2  +B2 -2BRsin(  9-£)-R 


(9) 


With  equations  (3)  to  (9)  the  topographic  height  z  at  each  point  on  the  y-axes  can  be  determined  as  well  as  all 
the  interrelations  between  all  parameters  for  one  Pixel. 


The  Interferogramme  for  the  whole  scene  is  generated  by  multiplying  pixel  by  pixel  the  image  7*i  with  the 
complex  conjugated  image  7*2  of  the  other.  Unfortunately,  the  phase  is  represented  as  modulo  2tt  only,  the 
phase  is  wrapped.  Thus  a  “phase  unwrapping”  has  to  be  carried  out  to  obtain  the  absolute  phase,  which  fi¬ 
nally  can  be  transformed  to  the  desired  height  z  at  each  point  over  the  x-y-area  leading  to  a  three  dimensional 
terrain  elevation  model. 


For  simplification  it  will  be  assumed  B«R  which  is  realistic  for  the  most  airborne  and  all  space-borne  SAR 
systems  used  up  to  now.  Under  that  condition  results  from  a  Taylor  development  of  (9) 


0»m 


2n 


2R 


Bsin(&-$) 


(10) 


For  both  space-borne  and  the  most  airborne  SAR  the  first  term  can  be  neglected  also  because 

^—«B  sin(i9-£) 

IK 


For  the  most  operational  SAR  Systems  results: 


2n 

0  »  —m  -y-  Bsin( ■ 9  —  £)  = 


B 

par 


(11) 


In  (1 1)  the  factor  Bpar  =  B  sin( &  -  £)  is  the  Baseline  component  parallel  to  the  look  direction. 

2.2  Interferometric  Performance,  Sensitivities  and  Errors. 

All  equations  content  both  Radar  System  parameter  and  Interferometry  parameter.  For  basic  considerations 
the  following  Radar  System  parameter,  principally,  are  assumed  to  be  known  a  priori:  The  height  H  above  a 
reference  area  (a  special  Geoid  ore  a  flat  area  respectively),  the  baseline  B,  the  wavelength  X  used,  the  plat¬ 
form  attitude  represented  by  the  baseline  tilt  angle  4,  the  look  angle  9,  and  the  slant  range  distance  R  be¬ 
tween  the  pixel  observed  at  P0,  the  both  antennas  A  /  and  A2  and  the  systems  Phase  Noise. 

However,  in  reality  all  these  parameters  have  to  be  determined  by  measurements  during  the  flight  and,  there¬ 
fore,  errors  will  be  introduced  due  to  unavoidable  measurement  inaccuracies.  That  will  be  demonstrated  ex¬ 
emplary  for  the  dependency  of  the  elevation  variation  Az  from  the  height  variation  AH.  Principally  holds 
under  the  assumption  of  small  pixels  for  the  variation  of  z  as  a  function  of  the  variation  AH  of  H: 

Az  dz  dz 

~  dH  ^  ^  *  m  AH 
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The  exemplary  starting  equation  is  (7)  which  gives  the  elevation  z  of  the  point  P  in  Fig.2.  In  connection  with 
(7)  results  for  the  elevation  sensitivity  Az: 

Az  =  AH  (13) 


From  the  sensitivity  Az  the  quadratic  error  az  results  in  that  case: 

(  dz 


< 7 \  = 


v  3!  I  j 


o\ 


The  entire  elevation  measurement  accuracy  depends  on  the  accuracy  of  the  determination  of  both  all  the  Ra¬ 
dar  System  parameters  and  all  the  Interferometry  parameters. 

The  general  expression  for  the  elevation  errors  is  derived  by  adding  the  all  quadratic  errors: 
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(14) 


The  subscripts  in  (14)  determine  the  respective  Interferometry  parameter. 


Special  interest  has  to  be  addressed  to  the  5th  term,  the  Phase  error  O®.  The  total  phase  results  from  the  sum 

(15) 


of  several  contributions.  The  designations  are  self  explaining: 

^flateartli^^ topography^^ 'atmosphere^® noise  Decorrelation 


The  Phase  errors  O motion  caused  by  uncontrolled  platform  motions  during  the  flight  will  be  discussed  more 
detailed  in  the  Systems  Lecture. 

The  sensitivity  of  a  function  0  =  0  (z,y)  against  variations  Az  and  Ay  respectively  results  from  differentia¬ 
tion  of  the  appropriate  equations  and  can  be  expressed  by  the  total  Differential 


dO  dO 

dO(z,y)=—dz+—dy 
dz  dy 


By  replacing  the  symbol  d  with  A  and  using  the  equations  (7),  (9),  and  (12)  results: 


(16) 


dO  dOd3  , dO  dOd3 

and 


^  do  .  do  . 

AO- - Az+ - Ay 

dz  dy  with  dz  d&  dz  dy  dS  dy 


From  (11)  follows: 


-~m^n  Bcos(  Q-E,)--m^n  B 
d3  A  A  eJJ 


The  factor  B  cos(T-q)  is  called  the  effective  Baseline  Beff  vertical  to  the  looking  direction. 


With  (6)  follows  both: 


Beff  =Bc°s(  3-%) 


d3  1  d3  1 

— = - and  — = - 

dz  R  sin  3  dy  R  cos  3 


(17) 


(18a)  and  (18b) 
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Therefore,  for  the  total  phase  variation  in  z-  and  y-direction  results: 

2k  B  cos(9 -£,)  2n  B  cos(9-%) 

AO  =  -m  — d  - ■ — n —  Az  -  m  —5 — ^  „ —  Ay  with  -  n  <  O  <  +k  mod  2n 

A  K  sind  A  K  cos  9  ■  (19) 

This  equation,  principally,  answers  the  question:  Which  Phase  variation  AO  will  be  caused  through  terrain 
variations  Az  in  vertical  ore  Ay  in  horizontal  direction?  It  delivers  the  interconnection  between  the  measured 
phases  and  the  height  characteristics  of  the  observed  pixel.  The  first  term  represents  the  phase  sensitivity 
against  elevation  variations  in  z-direction,  the  second  term  represents  the  sensitivity  against  a  horizontal  shift 
in  a  constant  altitude  z  in  y-direction  taking  into  account  the  relations  (3)  and  (6)  again.  Both  terms  represent 
interference  patterns. 

The  second  term  is  caused  from  the  change  in  slant  range  observation  distance  R  over  the  whole  swath  caus¬ 
ing  a  fringe  pattern  parallel  to  the  flight  direction  also.  Equation  (19)  shows:  the  flat  earth  alone  always 
causes  a  fringe  pattern  which  overlays  the  topographic  fringes  additionally  and,  in  order  to  observe  and  to 
analyze  the  topographic  fringes  only  these  two  effects  have  to  be  separated  by  subtraction  of  the  flat  earth 
interferogram  from  the  total  observed  interferometric  phase  (19). 

The  flat  earth  pattern  is  always  present  within  the  measurements  because  both  R  and  (R  +  r)  increase  with 
increasing  incidence  angle,  and  Beff,  principally,  changes  its  value  also.  The  flat  earth  pattern  has  to  be  re¬ 
moved  in  order  to  obtain  the  pure  elevation  above  the  reference  area  which  in  this  case  is  the  flat  area  at  z=0. 


AO 

z 


2n  Beff 

A  R  sin  9  ~  m 


2  k 

AT 


AS 


vert 


(20) 


2  k  Beff  2  k 

flat  A  R  cos  9  -  A  eff  nor  ^] ) 

The  relation  (20)  delivers  the  required  interference  pattern  between  the  measured  phases  and  the  terrain  ele¬ 
vation  z,  i.e.  a  vertical  shift  of  the  observed  point  P,  the  relation  (21)  delivers  the  interference  pattern  which 
results  from  a  horizontal  shift  of  the  observed  pixel.  AOvert  and  Aflhor  are  the  small  looking  angle  variations 
for  the  horizontal  and  the  vertical  shifts  of  the  point  P  respectively.  From  the  geometry  depicted  in  Fig.3  the 
following  relations  for  Aflverl  and  Aflh„r  can  be  determined: 

For  the  horizontal  shift: 

RA9  ,  »  R  sin(  AS  ,  ) 

nor  '  nor 


AR 

tan  9  (22) 
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For  the  vertical  shift: 


Az 

RA3  ,  »  R  sin(  A3  .  )  =  . - q 

vert  vert 7  tan  3 

Flere  with  result  from  (19)  and  (20)  the  phase  sensitivities  for  elevation  and  slant  range  variations. 

2n  Beff 

AO  r,  <  =  —m  -  n  ,  o  AR 
flat  A  R  tan  3 

Beff 

AO  =  —m  —x — „ — t — n  Az 

z  A  R  sin  3  _u<  (h  <tt 


(23) 


(24) 


(25) 


The  equations  (21)  and  (24)  represent  the  Interferogramme  of  the  flat  earth,  i.e.  an  Interferogramme  in  ab¬ 
sence  of  any  topographic  changes.  The  equations  (20)  and  (25)  represent  the  so  called  “flattened  topographic 
Interferogramme”  of  the  observed  area  in  slant  range  representation.  Both  Interferogrammes  have  to  be 
transformed  into  ground  range  presentations  in  order  to  obtain  the  desired  results.  The  interconnection  be¬ 
tween  ground  range  Rg  and  slant  range  R  is: 

Rg  =  R  sin  9 


With  that  relation  results  for  the  both  Interferogrammes: 


.  ,  2  it  AR 

AO  ,,  =-m — B  - 

1  eff  R 


cosi9 


g 


(26) 
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.  ,  2n  Az 

AO  —  —  m — B  re - 

z  X,  eff  r 

g 


(27) 


In  all  equations  varies  <D  mod  2n  between  -n  and  +n,  whilst  A<J)  varies  between  -n  and  +n  only.  The  width 
of  the  fringes,  principally,  is  given  for  AO  =  2n.  The  sensitivity  depends  on  the  minimum  phase  measure¬ 
ment  capability  AOmin  of  the  sensor.  The  state  of  the  art  is  AO>min  ~  5°  which  roughly  corresponds  to  1/70  of 
a  wavelength.  Fig.4  shows  simply  modeled  fringe  patterns  of  a  cone;  in  Fig.  5  a  real  measurement  result 
obtained  from  a  hill  formed  like  a  cone  is  depicted,  in  the  Figures  6,  and  7  for  a  space-borne  and  an  airborne 
case  measured  Interferogrammes  over  partially  relative  flat  areas  are  depicted. 


Fig.4  Models  of  fringe  patterns  for  a  flat  area  (left)  Fig.5  Measurement  results  over  a  flat  area  with  a 
and  a  modelled  cone  (right)  refuse  Cone,  amplitude  (top  left),  phase  fringes  (top 

right),  DEM  below.  The  shadowing  effect  in  the 
SAR  image  causing  a  complete  loss  of  coherence  is 
evident  in  the  amplitude  image 
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Fig.6  Set  of  phase  fringes  without  flat  earth  removal,  taken  simultaneously  in  X-,  C-,  and  L-Band  with  a  two 
pass  Interferometry  experiment  with  SIR-C  over  the  Mount  Etna  Test  area.  Time  difference  between  the  two 
passes  was  one  day.  The  Interferogramme  after  flat  earth  correction  is  depicted  in  Fig  20.  (page  21).  The  area 
in  the  left  part  was  relatively  flat,  the  largest  slopes  of  the  volcano  is  nearly  in  the  middle.  The  image  size  is 
66.8  km  in  azimuth  and  9  km  in  ground  range  (16000  radar  pulses  and  1024  range  bins  in  X-band  and  2048 
range  bins  in  L-  and  C-band).  More  results  of  the  same  scene  in  the  Fig.  18  to  24. 


Near  Range 

Far  Range 

$  =  25° 

o 

(ft 

u> 

II 

ef> 

R  =  3530  m 

R  =5579  m 

B  =  28.1  m 

11,5  m 

2 7i  -  10,3m 
Height 

2 n  -  49,6m 
Height 

2x  -  26,9  m 
Ground  Range 

2n  - 105,7  m 
Ground  Range 

Fig.7  Example  of  an  airborne  Interferogramme,  taken  with  E-SAR  in  L-band  repeat-pass  mode  in  the  rela¬ 
tively  flat  area  around  Oberpfaffenho fen/Germany.  Mainly  flat-earth  fringes  are  visible,  their  frequency 
changes  strongly  from  near-range  to  far  range. 
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For  the  space-borne  case  with  its  large  altitudes  and  small  incidence  angle  variations  a  nearly  periodical 
fringe  pattern  results  with  nearly  equal  fringe  width’s  result.  The  flat-earth  phase  is  a  linear  phase  ramp.  For 
airborne  data  the  situation  is  more  complicated,  the  look-angle  changes  strongly  from  near-range  to  far- 
range,  and  with  it  both  the  observation  distance  r  and  the  effective  base  line.  Therefore  nearly  everything  gets 
range-dependent.  The  flat-earth  phase  is  not  a  linear  phase  ramp  anymore,  it  has  a  curved  shape.  The  height 
resolution  changes  with  the  incidence  angle  also. 

2.3  The  Spectral  Approach  to  Interferometry 

The  relation  between  the  phase  of  SAR  Interferograms  and  ground  elevation  is  explained  as  usual  by  means 
of  geometrical  approaches  based  on  the  assumption  that  the  RF  bandwidth  is  so  small  to  be  negligible.  Thus, 
the  system  is  considered  monochromatic.  Flowever,  if  an  extended  system  bandwidth  W  is  introduced,  a 
relative  shift  of  the  ground  wave  number  spectra  dependent  on  the  baseline  and  the  local  slope  is  recognized. 
Important  consequences  are  revealed  from  this  result.  [Gatelli  et  al.  1 994] 

1)  It  can  be  shown  that  as  the  baseline  increases  and  in  the  absence  of  volumetric  effects  the  spatial  resolu¬ 
tion  of  the  Interferogramme  is  reduced  only  since  the  “geometrical  de-correlation”  can  be  removed  [Prati  C., 
Rocca  F.,  1993].  Therefore,  techniques  for  improving  the  Interferogramme  quality  (i.e.,  yielding  a  higher 
SNR)  can  be  found  (with  a  considerable  benefit  for  phase  unwrapping).  Moreover,  this  effect  could  be  com¬ 
pletely  avoided  by  means  of  tunable  systems. 

2)  Rules  for  the  design  of  pre-summing  filters  for  quick-look  interferometric  processing  can  be  derived. 

3)  Repeated  satellite  surveys  can  be  combined  in  order  to  improve  the  slant  range  image  resolution  because 
the  spectra  of  different  signals  correspond  to  different  bands  of  the  ground  reflectivity’s  spectrum.  The  rela¬ 
tion  between  the  frequency / and  the  ground  range  wave  number  ky  is  given  by: 

4  71  4  Tcf 

k  =  —  sin(9-a)=  -sin(9-a) 

y  A  c  (28) 


ORGANIZATION 


In  (28),  for  simplicity,  a  constant  uniform  slope  a  of  the  terrain  is  considered.  For  the  variation  of  ky  gener¬ 
ated  by  a  slight  change  of  the  looking-angle  Afi  follows: 


47JfA&  cos(S-a) 


(29) 


Thus,  in  general,  a  looking-angle  difference  Afi  generates  a  shift  and  a  stretch  of  the  imaged  terrain  spectra. 
However,  if  the  relative  system  bandwidth  is  small,  the  frequency /in  the  second  term  of  (28)  can  be  substi¬ 
tuted  with  the  central  frequency  f0. 


4nfQA9 

c 


cos(  9-a  ) 


(30) 


From  the  fact  that  the  radar  is  not  monochromatic,  it  has  a  bandwidth  W  centered  around  the  central  fre¬ 
quency  /o ,  results  that  by  changing  the  looking  angle  of  the  SAR  survey,  a  different  band  of  the  ground  re¬ 
flectivity  spectrum  shows  up. 


In  order  to  compare  the  shift  of  the  ground  reflectivity  spectrum  to  the  SAR  bandwidth  W,  the  ground  wave 
number  shift  of  (30)  is  expressed  as  an  equivalent  frequency  shift  A/  The  following  expression  of  A/  for  an 
angular  separation  All  can  be  obtained  by  differentiating  (28)  directly  and  taking  into  account  the  relations 
for  the  small  incidence  angle  changes  resulting  from  Fig.  3: 


Af=~ 


foA 9 


cB 


eff 


tan(  9-a  )  Rltan(  3-a  ) 


(31) 


Note:  (31)  does  not  state  that  by  changing  the  looking  angle  of  the  SAR  survey  the  radar  bandwidth  is  shifted 
by  A /  it  just  says  that  by  changing  the  SAR  looking  angle,  the  backscattered  signal  contains  different  spec¬ 
tral  components  of  the  ground  reflectivity  spectrum.  The  signals  received  by  two  SAR’s  separated  by  an 
angle  A9  have  the  same  spectral  components  of  the  first  signal  whilst  in  the  second  signal  the  whole  spec¬ 
trum  is  shifted  only  by  A f 
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A  similar  result  can  be  demonstrated  for  the  bistatic  SAR  case  where  a  stereo  pair  is  acquired  simultaneously 
by  two  sensors:  the  first  is  equipped  with  a  transmitter  and  a  receiver,  whereas  the  second  is  equipped  with  a 
receiver  only.  It  can  be  shown  that  the  spectral  shift  induced  by  a  baseline  Beff  is  half  that  obtained  in  the 
monostatic  case: 


4f=- 


2RAtan( 9-a  ) 


(32) 


In  the  mono  static  case  the  two  SAR  spectra  become  totally  disjoint  as  the  frequency  shift  A/equals  the  sys¬ 
tem  bandwidth  W.  The  critical  baseline  can  be  computed  as: 


B 


effcrit 


WRAtan(  9-a) 
c 


(33) 


By  using  the  geometric  range  resolution  5rge  instead  of  the  Bandwidth  W  results  from  the  well  known 
interconnection  between  range  resolution  and  bandwidth: 


B 


effcrit 


RXtan(  9 -a) 


(34) 


The  Tables  1  to  3  summarize  the  results  obtained  in  this  Chapter,  Table  4  is  an  example  for  real  results. 


Tab.  1  Basic  Interferometry  Formulas 


Ground  Range: 

y  =R  sinfi 

(6) 

Local  Terrain  Height: 

r2-B2 

*7  U 

„  Q 

2Bsin(9-^)-2r 

(7) 

Critical  Baseline 

WR2tan(  9-a  ) 

B  effcrit  c 

RAtan(  9-a  ) 

28 

rge 

(34) 

Surface  Change  Detection 
Differential  Interferometry 

X 

Az  =  j—  AO 

4  n 

Tab.  2  Parameter  Examples  for  Interferometry  with  a  space-borne  SAR  from  ERS/ASAR  Type 


Input  parameters 

Value 

•  velocity:  v 

7.464  km/s 

•  range  delay  (mid  swath):  T0 

5.779  ms 

•  slant  range  (mid  swath):  R0 

866.656  km 

•  incidence  angle  (mid  swath):  9 

23° 

•  altitude:  H 

794  km 

•  wavelength:  X 

0.0566  m 

•  baseline  tilt  angle: 

0° 

•  baseline  length:  B 

200  m 

•  surface  slope:  rx 

rv 

•  phase  error:  AO 

5° 

•  baseline  length  error:  AB 

0.001  m 

•  baseline  tilt  angle  error:  Aq 

1° 

•  slant  range  error:  A R 

3  m 

•  height  error:  AH 

1  m 

RTO-EN-SET-O81bis 


2  - 11 


Normal  and  Differential  SAR  Interferometry 


ORGANIZATION 


Tab.  3  Resulting  Errors  for  absolute  Height  Estimation 


Error 

Type 

Formula 

Errors 

Phase 

2  Tf  •  Q 

AR  sin  .9  g 

^ 0  ~  4xBcos(S-%)  A0  ~  4nBeff  A° 

Ad>=5°  ^  Az®  0,723m 

Baseline 

Length 

AB 

Az  g  =  -R  tan(  &  -  £,)  sin  <9  -g- 

AB  =  0,001  m  ^,Azb=  0,719  m 

Baseline 

Tilt  Angle 

Az  ^  =  ~Rg  sin 

A^  =  1°  Az^  =5910,20  m 

Range 

Delay 

Az  =  -  cos  9AR 

A R  =  3m4  AT  =  0,2  psec 
^  Azr=  -2,762  m 

Altitude 

Azr  =  AH 

A H  =lm4  Azh=  lm 

Tab.  4  Digital  Elevation  Model  errors  caused  by  instrument  inaccuracies  of  SRTM/X-SAR.  Goal  was 
a  relative  height  accuracy  (90  %)  <  6  m  absolute  height  accuracy  (90  %)  <  16  m  [Keydel  et  al  2000] 


DEM  Error  Examples  for  SRTM  (Middle  of  Swath) 

Relative  (30  seconds) 

Absolute  (11-days) 

Error  Type 

Inaccuracy 

DEM-Error 

Inaccuracy 

DEM-Error 

Baseline  Tilt  Angle 

2  arcsec 

3,0  m 

9  arcsec 

13,4  m 

Baseline  Length 

1,3  mm 

0,8  m 

4,0  mm 

2,6  m 

Instrument  Phase 

4,0deg 

4,2  m 

4,0  deg 

4,2  m 

Total  (RSS) 

5.5  m 

14.4  m 

Conclusions 

Absolute  height  measurement  with  interferometric  SAR  requires  an  extremely  stable  system  performance. 
Especially  the  orbit  has  to  be  known  very  precise  due  to  its  great  impact  on  the  baseline  length,  the  tilt  angle 
accuracy  and  the  altitude  accuracy.  An  instrument  phase  error  of  5°  leads  to  an  error  in  the  height  estimation 
of  about  1.5  m  in  the  mid  swath  and,  therefore,  seems  to  be  acceptable. 


3.  Coherence  and  Correlation  in  SAR  Interferometry 

For  SAR  Interferometry  the  coherence  between  two  signals  to  be  compared  defined  in  the  Introduction  Lec¬ 
ture  is  a  fundamental  parameter.  The  Coherence  is  a  measure  for  the  quality  of  Interferogrammes.  Equation 
(4)  will  be  used  to  produce  the  complex  SAR-Interferogramme  replacing  the  Fields  Ei  and  E2by  the  complex 
Images  E  and  E  gained  in  each  of  the  two  Interferometry  channels.  This  Interferogramme  is  defined  as  the 
product  of  the  complex  SAR  values  of  the  second  image  with  the  complex  conjugate  of  the  reference  image. 
Principally,  the  Interferogramme  is  the  complex  Coherence  Image  Cimage: 

r  {(/,(»/;((») 

(<))(/,  (o/:<o) 
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The  corresponding  amplitudes  have  to  be  averaged  and  the  corresponding  phases  have  to  be  differenced  at 
each  point  in  the  image.  The  resultant  phase  of  the  Interferogramme,  which  is  directly  related  to  the  topogra¬ 
phy,  is  measured  modulo  2n  only.  To  calculate  the  elevation  of  each  point  it  is,  therefore,  necessary  to  add 
the  correct  integer  number  of  phase  cycles  to  each  phase  measurement.  The  method  of  solving  this  ambiguity 
is  called  phase  unwrapping. 

On  that  basis  applications  of  interferometry  can  be  split  into  two  categories,  based  on  whether  the  main  pa¬ 
rameter  of  interest  in  the  complex  interferograme  is  the  phase  (4)  or  the  magnitude  (5)  of  the  complex  coher¬ 
ence,  Fig  8. 


Fig.8  SIRC  result,  L-Band.  Example  for  fringe  degradation  due  to  coherence  loss  caused  by  temporal  de- 
correlation  in  the  Mount  Etna  area  which  is  covered  by  vegetation  in  the  yellow  Square. 

The  phase  information  can  be  used  for  Digital  Elevation  Model  (DEM)  production,  provided  the  scene  is 
relatively  stable.  For  repeat-pass  Interferometry,  the  phase  may  also  be  used  for  studying  small  coherent 
motions  of  the  target  between  the  imaging  times,  a  motion  in  which  most  of  the  scatterers  move  simultane¬ 
ously.  This  finds  applications  in  geological  studies  of  land  motion  (e.g.  associated  with  seismic  activity  or 
land  slip),  and  in  glaciology.  The  accuracy  of  the  measurements  of  height  and  motion  are  in  part  determined 
by  the  noise  of  the  phase  measurements.  This  phase  noise  is  usually  characterized  by  the  magnitude  of  the 
coherence  of  the  Interferogramme,  often  simply  referred  to  as  the  coherence. 

Important  information  can  be  extracted  from  images  using  the  coherence.  This  is  true  for  land  use  classifica¬ 
tion,  where  the  information  obtained  is  complementary  to  the  intensity  of  the  SAR  images,  and  may  give,  for 
example,  better  discrimination  between  clear  cuts  and  forest  areas  than  intensity  measurements.  Another  use 
of  coherence  in  this  respect  is  the  monitoring  of  changes  to  land  surfaces  (e.g.  harvesting  of  agricultural 
crops).  Fig  9  shows  an  example  for  Land  classification  obtained  with  coherence  evaluation  of  a  multi  fre¬ 
quency  X-SAR/SIRC  result;  the  respective  phase  image  of  that  very  scene  without  flat  earth  removal  is  de¬ 
picted  in  Fig.6  (page  8). 
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Coherency  Classification  Mop 
based  on  X-  C—  and  L-bond 

Low  and  spare*  vegetation 

Fresh  osh  and  scoria 

Otd  osh  ond  scoria  mantle  A 

Ok)  covered  kjva  flows 

Clear  lava  flows 


Test  site:  Mt.  Etno  /  Sicily 


Fig.  9  SIR-C  Classification  result  of  the  Mount  Etna  test  area  obtained  with  coherence  comparisons  in  L-,  C-, 
and  X-Band 

There  are  a  number  of  causes  for  de-correlation  in  Interferometry,  and  it  is  possible  to  write  the  coherence  as 
a  product  of  factors  which  give  some  physical  insight  into  the  origin  of  the  de-correlation: 


Y  Yinstrument  YprocessorYbaselineYazimuthYnoise  Ytemporal  Yspatial 

The  first  three  terms  are  totally  independent  of  the  target  being  imaged,  the  fourth  and  fifth  are  related  to 
both  the  system  and  the  target,  and  the  final  term  is  totally  dependent  on  the  target  properties  as  will  be  ex¬ 
plained. 

Yinstrument  refers  to  both  Radar  and  Platform 

Yprocessor  refers  to  effects  produced  during  the  processing  of  the  SAR  images  to  form  the  Interferogramme 
and  includes,  for  example,  miss-registration  of  the  images  when  aligning  them  before  forming  the  Interfero¬ 
gramme  and  overlap  of  the  range  and  azimuth  spectra.  The  latter  means  the  two  surveys  must  image  the 
scene  under  the  same  squint  angle,  they  have  to  have  the  same  Doppler  Centroid.  The  to  days  processing 
routines  are  estimated  to  register  the  two  images  to  within  one  hundredth  of  a  pixel  in  both  range  and  azi¬ 
muth,  leading  to  a  negligible  loss  of  coherence. 

Yazimuth  is  a  factor  introduced  to  account  for  differences  in  the  Doppler  Centroid  frequencies  between  the 
two  images.  For  repeat-pass  Interferometry  exemplary  with  the  ERS-1  or  RadarSAT-1  satellites  this  effect  is 
negligible  as  the  high  stability  means  that  the  Doppler  centroids  typically  vary  by  only  about  0.02  of  the 
sampling  frequency. 

Ynoise  represents  the  effect  of  noise  in  the  system  is  well  known  and  can  be  expressed  as 

CNR  1 

y"oise  =  CNR+1  =  1+CNR 

where  CNR  denotes  the  clutter-to-noise  ratio.  This  is  important  for  regions  of  weak  backscattering. 

Ytemporal  refers  to  the  de-correlation  in  single  pass  Interferometry  which  depends  simply  on  the  factors  dis¬ 
cussed  above.  In  the  case  of  repeat  pass  Interferometry,  where  the  time  interval  between  images  may  be  days 
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or  longer,  there  is  also  a  contribution  to  the  de-correlation  due  to  the  possibility  of  changes  in  the  scatterers 
themselves  and  due  to  the  changes  of  weather  conditions  in  the  propagation  path.  An  agricultural  field  for 
example  de-correlates  steadily  with  the  evolution  of  the  crop.  The  correlation  is  completely  destroyed  once 
the  field  is  ploughed.  Forests  de-correlate  even  more  rapidly  etc.  On  the  other  hand  many  manmade  and 
natural  objects  prove  to  stay  coherent  over  years  and  decades  which  enables  the  techniques  of  Permanent 
Scattering  Interferometry  considered  in  Chap. 7 

Ybaseline  refers  to  the  baseline  of  the  two  passes  which  has  to  be  smaller  than  the  critical  effective  baseline 
Yspatial=  Yrange  Yvolume  refers  to  the  de-correlation  caused  by  different  acquisition  look  angles 
Yrange  refers  to  the  de-correlation  due  to  different  projections  of  the  bandwidth  onto  ground 
Yvolume  refers  to  the  de-correlation  related  to  the  height  distribution  of  the  scatterers 

In  Fig.  8  the  influence  of  coherence  degradations  on  the  Interferogramme  is  depicted.  It  is  impossible  to  ap¬ 
ply  the  definition  of  coherence  directly  to  the  SAR  images,  due  to  the  fact  that  it  involves  ensemble  averag¬ 
ing.  In  practice  an  assumption  is  made  in  which  a  spatial  average  is  taken  to  be  a  reasonable  approximation 
of  an  ensemble  average.  This  requires  that  phase  variations  between  the  pixels  being  averaged  have  been 
corrected  for  topographic  variations.  The  averaging  process  introduces  a  bias  into  the  coherence  measure¬ 
ments,  the  level  of  which  is  determined  by  the  number  of  independent  samples  used. 

However,  the  coherence  of  a  pair  of  interferometric  images  gives  one  indication  of  the  quality  of  the  ob¬ 
tained  Interferogramme  [Monti  97].  Also  the  coherence  matrix  is  a  basic  starting  point  for  some  unwrapping 
algorithms.  If  real  time  is  desired,  or  an  increment  in  the  speed  of  all  the  chain  from  raw  data  to  the  digital 
elevation  model  (DEM)  is  needed,  the  time  reduction  in  computing  the  coherence  should  be  undertaken. 

The  usual  coherence  estimator  needs  to  calculate  the  local  frequency  of  the  interferometric  fringes  very  accu¬ 
rately  and  therefore  this  algorithm  becomes  very  time  consuming.  The  coherence  between  two  complex  SAR 
images  I ,  and  I2  is  given  by  (28).  Assuming  stationary  regions  and  ergodicity  (i.e.  all  scatterers  within  the 
estimation  area  are  independent  with  identical  statistical  properties),  the  coherence  can  be  calculated  with 
high  accuracy  as: 


r 


est 


N  M 

Z  Z  1 1/  \e 

n  n  J(n>m)  2(n,m ) 
n=0  m=0 


-j@(  n,m) 


N  M 

Z  Z 

n = 0  m = 0 


I j  (  n,m) 


2  N  M 

Z  z \l2(n,m) 
n=0 m=0 


(29) 


yest  is  called  coherence  estimator.  N  and  M  give  the  size  of  the  estimation  window  in  range  and  azimuth.  The 
factor  The  inclusion  of  a  factor  of  is  due  to  the  difference  between  the  two  SAR  images,  which  differ 

by  the  deterministic  phase,  that  should  be  estimated  and  corrected. 


4.  Generation  of  Interferometric  Products 
4.1  Interferometric  Product  Generation  Principle 

The  input  data  used  for  SAR  Interferometry  are  the  complex  SAR  image  pairs  and  the  respective  navigation 
data.  The  procedure  used  for  the  generation  of  the  interferometric  products  is  depicted  in  Fig.  10  and  can  be 
summarized  as  follows: 
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Fig.  10  Principal  procedure  used  for  the  generation  of  interferometric  products,  exemplary  DEM’s  and  co¬ 
herence  maps. 

1.  Processing  including  range  and  azimuth  fdtering  of  the  complex  images.  Before  the  Interferogramme  is 
formed,  the  complex  image  pair  has  to  be  filtered  separately  in  both  range  and  azimuth  directions  to  optimize 
the  coherence  and  the  interferometric  phase  purity.  The  so-called  spectral  shift  range  fdtering  will  be  applied 
in  order  to  minimize  baseline-induced  de-correlation.  Filtering  to  a  common  Doppler  band  in  the  azimuth 
direction  will  be  applied  to  eliminate  the  de-correlation  effect  induced  by  the  Doppler  Centroid  difference 
between  the  two  images  [J.  Moreira,  M.  Schwab isch,et.al  1995]. 

2.  Coregistration.  In  order  to  obtain  the  interferometric  phase,  the  image  pair  has  to  be  co-registered.  This 
has  to  be  carried  out  with  a  high  degree  of  accuracy;  state  of  the  art  is  the  order  of  one  hundredth’s  of  a  pixel. 
3a.  Interferogramme  generation.  The  Interferogramme  is  generated  by  multiplying  the  first  image,  L  ,  with 
the  complex  conjugate  of  the  second  image,  referred  to  as  h 

3b.  Flat-earth  phase  removal.  To  perform  the  Interferogramme  noise  filtering  and  simplify  the  phase  un¬ 
wrapping  procedure,  a  first-order,  flat-earth  phase  removal  is  carried  out.  The  flat-earth  phase  is  computed 
by  measuring  the  dominant  fringe  frequency  in  the  azimuth  and  range  directions  of  the  Interferogramme.  It 
thus  consists  of  two  different  linear  phases  in  the  range  and  azimuth  directions.. 

3c.  Phase  noise  reduction.  In  order  to  reduce  the  phase  noise  and  to  facilitate  the  phase  unwrapping,  the 
Interferogramme  has  to  be  filtered  using  an  averaging  window  with  a  size  of  several  resolution  elements  in 
both  range  and  azimuth.  Respective  pixel  spacings  and  resolutions  used  for  X-SAR  processing  are  shown  in 
Table  1  exemplary.  After  the  filtering,  the  interferometric  phase  is  determined  by  calculating  the  arctangent 
of  the  Interferogramme.  The  terminology  and  respective  procedures  are  given  in  the  next  sections. 

4.  Phase  Unwrapping,  the  most  important  procedure  will  be  deeper  described  in  subchapter  4.2 

5.  Post  Processing  with  height  and  velocity  estimation  ending  up  with  the  desired  DEM  and  a  coherency 
map  of  the  scene. 


4.2  Phase  Unwrapping 

Phase  Unwrapping  is  the  key  procedure  for  SAR  Interferometry  data  evaluation  with  respect  to  establishing 
DEM’s.  Meanwhile,  there  are  many  methods  for  the  two  dimensional  Phase  unwrapping.  The  very  basic 
principle  is  depicted  in  Fig.  1 1  for  a  simple  one  dimensional  case  which  is  self  explaining. 
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Fig  1 1  Basic  principle  of  the  phase  unwrapping  procedure,  Reigber  2001 

A  more  sophisticated  principle  is  the  least  square-Method  which  minimizes  the  differences  between  the  de¬ 
rivatives  of  the  wrapped  and  unwrapped  phase  functions  using  the  fact  that  the  minimum  of  the  first  deriva¬ 
tive  is  determined  by  the  null  of  the  second  derivative.  In  the  two  dimensional  case  the  partial  differential 
quotients  have  to  be  used.  This  method  seems  to  be  very  simple  and  fast.  However,  the  errors  are  propagat¬ 
ing  very  strongly  down  to  the  entire  sequence  and  therefore,  this  method  is  not  very  practicable  for  large  area 
processing. 

As  further  examples  branch  cutting,  fringe  detection,  and  knowledge  injection  should  be  mentioned  in  a  bit 
more  detail  despite  of  the  fact  that  they  don’t  represent  the  most  modem  state  of  the  art.  However,  they  seem 
to  be  important  to  gain  insight  into  the  whole  complex.  [Gens  R.,  1996] 

Branch  cutting 

The  first  method,  which  is  based  on  creating  so  called  branch  cuts  was  developed  in  1988  already  by  Gold¬ 
stein  el  al.  They  pointed  out  that  there  are  two  possible  types  of  errors  in  the  unwrapped  sequence.  The  first 
are  local  errors  in  which  only  a  few  points  are  corrupted  by  noise;  the  second  are  global  errors  in  which  the 
local  error  may  be  propagated  down  the  entire  sequence.  The  global  errors  arise  from  the  residues.  Residues 
are  local  errors  in  the  measured  phase  caused  by  the  noise  in  the  signals  or  by  actual  discontinuities,  i.e., 
layovers,  in  the  data.  There  are  three  major  sources  of  noise.  The  first  is  thermal  noise,  the  second  is  the 
speckle,  which  is  a  random  variable,  and  the  third  one  is  caused  by  an  antenna  effect  depending  on  the  an¬ 
tenna  separation,  the  bandwidth  of  the  radar,  and  the  geometry  of  the  image  formation.  Assuming  that  the 
original  scene  is  sampled  often  enough  so  that  the  true  phase  will  not  change  by  as  much  as  one -half  cycle 
per  sample  point  the  phases  can  be  unwrapped  in  a  consistent  manner.  This  can  be  done  if  the  residues  are 
identified  and  suitable  branch  cuts  are  made  between  the  residues  to  prevent  any  integration  path  from  cross¬ 
ing  these  cuts.  Filtering  the  signal  will  reduce  the  residues  in  the  image  but  it  also  reduces  the  spatial  resolu¬ 
tion  in  the  final  topographic  map.  Before  integrating  the  phase  differences  all  residues  which  were  previously 
low-pass  filtered  are  to  connect  with  branch  cuts.  The  phase  differences  are  then  integrated  without  crossing 
any  of  the  cuts. 
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Fringe  detection 

Another  method  for  the  two-dimensional  phase  unwrapping  is  the  fringe  detection  method  of  Lin  et  al. 
(1992).  The  first  step  in  this  method  is  to  find  the  location  of  the  fringe  lines  in  the  phase  difference  image  by 
using  edge  detection  techniques  which  mostly  include  the  two  steps  of  enhancement  by  filtering  and  then 
applying  a  threshold  to  the  image.  The  phase  difference  image  is  median  filtered  with  a  window  size  which 
depends  on  the  density  of  the  fringes.  For  the  detection  of  the  fringe  lines,  five  edge  masks  with  different 
orientations  are  applied.  The  phase  is  then  unwrapped  by  adding  a  multiple  of  2n  each  time  the  integration 
path  crosses  a  fringe  line.  The  multiple  depends  on  how  many  fringe  lines  lie  between  the  pixel  and  a  refer¬ 
ence  pixel.  The  solution  is  considered  good  when  the  fringes  are  greatly  separated  and  the  signal-to-noise 
ratio  (SNR)  is  high  in  the  entire  wrapped  phase  image  [Flartl  et  al.l 993]. 

Knowledge  injection 

This  method  is  used  for  reducing  the  density  of  fringes  by  incoiporating  the  use  of  a  rough  DEM.  [Adragna 
1995].  The  same  technique  is  also  used  for  flat  Earth  removal,  and  for  differential  Interferometry  (Chap. .5.), 
in  which  knowledge  of  the  topography  is  injected  to  isolate  the  terrain  movement.  It  should  be  noted  that  all 
the  methods  of  phase  unwrapping  described  above  are  generally  performed  after  flat  earth  removal. 

4.3  Co-Registration 

The  co-registration  is  a  very  important  processing  step  in  Interferometry  also.  For  distributed  scatterers,  a  co¬ 
registration  inaccuracy  introduces  substantial  noise  into  the  interferometric  phase  (Bamler  et  al.,  1993),  and 
point  scatterers  have  to  deal  with  a  systematic  phase  error  (Flolzner  et  al.,2001).  (Compare  the  PSI  considera¬ 
tions  in  Chap. 7  also).  Bamler  (2000)  showed  that  the  error  of  the  co-registration  parameter  estimation  for 
distributed  scatterers  depends  on  the  number  of  correlation  points  n  and  the  coherence  y: 


a 

corr 


2 


The  dimension  here  is  the  unit  of  the  resolution  cell.  Accordingly,  in  order  to  guarantee  appropriate  co¬ 
registration  accuracy  in  low  coherence  scenes  the  correlation  window  size  needs  to  be  increased.  Flowever, 
the  co-registration  still  strongly  depends  on  the  temporal  coherence  and  the  co-registration  limit  is  often  in¬ 
sufficient. 


5.  Along-Track  Interferometry  for  Velocity  Determination 

The  arrangement  of  the  two  antennas  in  cross  velocity  direction  is  essential  for  topographic  mapping  because 
that  arrangement  only  leads  to  phase  differences  between  the  two  SAR  images.  An  along  track  arrangement 
as  depicted  in  Fig.2  (page  3)  would  produce  two  identical  images  with  no  phase  differences  at  all  due  to  the 
same  processing  procedure.  Each  object  had  both  the  same  distance  and  the  same  Doppler  Frequency  with 
respect  to  the  two  receiver  ore  antennas  respectively.  Flowever,  there  is  a  time  shift  T  between  the  two  im¬ 
ages  due  to  the  fact  that  both  antennas  pass  the  point  of  nearest  approach  to  each  pixel  -  which  is  the  locus  of 
Doppler  Frequency  0  -  in  a  sequence  which  is  dependent  on  Baseline  length  and  platform  velocity  vpiat. 

T=B/vp,at 

A  target  changing  its  position  during  that  time  with  the  velocity  vtarg  will  change  its  position  in  the  second 
image  in  accordance  with  its  own  velocity  producing  a  phase  difference  ®.  The  Phase  difference  between 
corresponding  pixels  is  a  measure  for  the  velocity  of  the  object  in  the  scene.  This  procedure  is  called  Along 
Track  Interferometry. 


2  n  2  7Z  t  arg 

0  =  m  v .  T  —  m  — s - B 

1  1  ”  plat 


(30) 


In  order  to  determine  the  ground  range  velocity  Vtargy  of  a  target  ore  a  pixel  respectively  which  is  moving  on 
ground  neglecting  the  vertical  velocity  component  vtarg  z  for  simplicity  reasons  results: 
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V 


targ 


targy  sin  9 


(31) 


The  measurement  of  the  velocities  of  different  resolution  cells  of  water  surfaces  for  instance  allows  the  esti¬ 
mation  of  water  currents  with  respect  to  both  current  velocity  and  the  respective  direction.  However,  the 
along-track  technique  will  not  be  suitable  for  the  measurement  of  along  track  velocity  components.  But  with 
two  separate  measurements  in  cross  directions  a  two  dimensional  current  map  can  be  established.  Fig.  12 
shows  an  example.  In  Fig.  13  results  of  moving  target  detection  and  measurement  results  are  depicted.  Here, 
additional  a  priori  map  information  with  respect  to  position  and  direction  of  the  respective  roads  within  the 
scene  has  been  used  assuming  that  the  identified  vehicles  were  driving  on  roads  only. 


Water  Current  Mapping: 

Interferometric  phase  corresponds  to  flow  velocities. 

to 

I 


Bottom  Image: 

Overlay  of  SAR  image  and  E-SAR 
interferometric  phase 

Colours  represent  flow  velocities: 
yellow:  2  m/s;  blue:  no  current 


Accuracy:  10  cm/s 


Application  Development  of  a  model 
to  predict 

sand  sedimentation 
based  on  water  movements 


0  m/sec 


2  m/sec 


Fig.  12  Example  for  current  velocity  measurements  with  along  track  interferometry.  The  colours  represent 
flow  velocities,  :yellow:  2  m/s;  blue:  no  current 
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Fig.  13  Examples  for  MTI  with  space  borne  (left)  and  airborne  (right)  along  track  Iinterferometry. 

Left:  X-SAR/SRTM  one  pass  Interferometry  image,  baseline  7m.  Motorway  A9,  Munich  -  Nuremberg, 
Germany,  motorway  blue,  red  arrows  vehicles  with  their  velocities  (yellow)  ,  black  arrow:  SRTM  track 
[Breit  et.al.2003].  Right:  airborne  one  pass  Interferometry  image  with  vehicles  in  two  directions  on  the  high¬ 
way  [Ender  1996] 


6.  Differential  SAR  Interferometry 

Many  operational  monitoring  tasks  require  the  detection  of  changes  during  a  certain  time  period.  This  can  be 
performed  with  Differential  SAR  Interferometry  which  registers  the  respective  topographic  changes  of  a 
DEM  which  happened  during  the  time  period  since  the  DEM  has  been  established.  There  are  principally  two 
techniques  depending  on  the  availability  of  pre-existing  digital  topographic  maps. 

If  such  a  DEM  with  acceptable  accuracy  is  available  then  it  has  to  be  subtracted  from  the  DEM  which  results 
from  the  interferometric  phase  field  measurements.  The  difference  between  the  two  DEM’s  represents  the 
changes  to  be  determined. 

If  no  digital  topographic  maps  are  available  then  two  Interferogrammes  of  the  same  area  have  to  be  taken 
over  the  desired  time  period.  The  first  one  is  used  to  generate  a  reference  Interferogramme  and  the  second 
one  to  generate  a  Interferogramme  which  contents  both  the  reference  DEM  together  wit  the  changes  to  be 
monitored.  This  leads  to  a  double-difference  interferogram,  which  can  be  derived  from  the  two  correspond¬ 
ing  Interferogrammes  by  subtraction  of  both  Interferogrammes. The  phase  discontinuities  in  the  two  Inter¬ 
ferogrammes  will  not  occur  in  the  same  places.  Therefore,  each  Interferogramme  must  be  unwrapped. 

The  sensitivity  of  SAR  Interferometry  depends  on  the  geometry  of  the  observation  and  can  be  increased  at 
the  expense  of  spatial  resolution  by  averaging  over  more  pixels.  The  length  of  the  baseline  in  relation  to  the 
wavelength  used  has  an  influence  on  the  amount  of  noise  which  is  introduced.  A  decreasing  baseline  ore  an 
increasing  wavelength  leads,  therefore,  to  a  higher  accuracy  of  the  detected  height  differences. 
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The  principal  geometry  for  the  space  borne  two  pass  mode  is  depicted  in  Fig  14.  For  simplicity  reasons  in 
describing  the  principles  it  is  assumed  that  both  Interferogrammes  use  the  same  primary  flight  path,  however 
they  can  have  different  baselines.  The  following  Interferogrammes  will  be  obtained: 


Fig.  14  Differential  Interferometry  Principle 


Interferogram  1: 


2n  2n 

0  r  =  m  rr  =  0  ,  ,  =  m  B , 

1  X  1  topography  A  1  par 


(32) 


Interferogram  2: 


2n 
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(33) 


The  subtraction  of  (32)  from  (33)  leads  to  the  displacement  in  line  of  sight  direction 


02  @2  @  displacement 


2tc 

=  m  X  P 


(34) 


This  expression  is  not  dependent  on  the  topography  anymore.  As  explained  in  Cap.  2.2  the  sensitivity  results 
from  differentiation  which  leads  to 


*0  2 
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Fig  15  shows  exemplary  the  displacements  caused  by  an  earthquake  in  Japan. 


Izmit/Turkey.  17  August.  1999 


Fig.  15  Differential  Interferometry  result.,  displacement  pattern  caused  by  an  earthquake 
[http://comet.nerc.ac.uk/schools_eq_izmit.html  and  Bamler,  2000] 


7.  Permanent  Scatterer  Interferometry  (PSI) 

Differential  Interferometric  SAR  over  longer  time  periods  is  often  limited  by  temporal  de-correlation.  As 
mentioned  already  many  manmade  and  natural  objects  prove  to  stay  coherent  over  years  and  decades.  These 
so-called  Permanent  Scatterers  (PS)  can  be  used  for  the  ultimate  sensitivity  in  motion  measurement  from 
space  (Ferretti  et  al.,  2000;  Ferretti  et  al.,  2001)  by  using  them  as  points  of  reference.  PS’s  are  usually  bright 
points  in  the  image,  with  a  very  high  signal-to-noise  ratio.  Flowever,  the  phase  information  itself  cannot  be 
used  for  detection  of  stable  points,  since  it  contains  phase  uncertainties  induced  by  unknown  elevation  and 
deformation  but  amplitude  time  series  can  be  well  used  for  this  [Ferretti  2001].  The  permanent  Scatterer 
technique  uses  a  sequence  of  many  Interferometric  SAR  products  of  the  same  area  of  interest  taken  over  a 
long  time  period  establishing  in  a  first  step  a  data  base  of  very  many  SAR  images.  From  the  observed  phase 
time  series  of  the  differential  Interferogrammes  the  elevation  difference  and  deformation  rate  difference  over 
time  can  be  determined. 

Presently  such  databases  are  available,  principally,  from  the  most  areas  of  the  globe  taken  from  different 
operational  satellites  like  ERS-1/2,  RadarSAT-1,  J-ERS  etc.  which  where  in  orbit  permanently  since  more 
than  15  years.  (Some  of  them  will  be  numbered  and  considered  in  the  Systems  chapter  of  this  very  lecture 
Series). 

If  several  SAR  images  are  supeiposed,  even  the  propagation  errors  due  to  tropospheric  water  vapour  can  be 
averaged  out.  In  practice  30  -  100  SAR  images  of  the  same  areas  are  used  to  form  an  interferometric  stack 
spanning  a  time  period  of  typically  5-12  years.  In  this  co-registered  stack,  the  PS’s  are  identified  and  their 
phase  history  can  be  tracked  and  analysed.  Several  filtering  and  adjustment  steps  are  required  to  remove  orbit 
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errors,  long-wavelength  ionospheric  delays  and  to  reduce  atmospheric  effects  and  noise.  The  result  is  a  map 
of  subsidence/inflation  rates  to  an  accuracy  of  less  than  1  mm/year.(Bamler,  SRTM  and  beyond) 

If,  for  example,  n  +  1  ERS  SAR  images  from  one  certain  area  are  available  the  data  will  be  first  co-registered 
on  a  unique  master  and  a  first  DEM  of  the  area  will  be  established  as  a  reference.  Then,  n  differential  Inter- 
ferogrammes  between  all  SAR  images  and  the  master  will  be  computed.  After  DEM  compensation,  the  re¬ 
sidual  phase  O/of  the  Interferogramme  number  i  is: 


0  . 
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noise  topo  displacement 

+  0  .  +  0  .  +  0  . 
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i  =  l..n 


(36) 


r,  (T,)  represents  the  momentary  target  position  at  the  time  T,  and  herewith,  principally,  the  possible  target 
motion  in  the  line-of-sight  direction,  T,"""  is  the  atmospheric  phase  contribution,  (fr  de"""  the  contribution  of 
the  de-correlation  noise,  <b  ‘opo  is  the  topographic  phase  error  etc.  All  errors  are  mentioned  already  in 
Chap.2.2, 


A  time  series  analysis  of  the  first  term  in  equation  1  reflects  the  target  motion:  here  for  simplicity  a  constant 
velocity  model  is  supposed.  An  uniform  strain  rate  hypothesis  is  often  used  in  geophysical  modelling;  how¬ 
ever,  more  complex  models  can  be  adopted  if  the  linear  model  is  unrealistic.  Permanent  scatterers  are  charac¬ 
terized  by  low  phase  residues  and  are  determined  by 
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The  first  term  in  equation  (36)  can  be  written  as  follows: 
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Vr  is  the  unknown  component  of  the  target  velocity  m  in  the  Line  of  Sight  (LOS)  direction  and  7j  is  the  tem¬ 
poral  baseline  between  the  master  acquisition  and  the  generic  i-th  slave  image.  Since  there  are  n  differential 
Interferogrammes  of  the  same  area  with  different  temporal  and  geometric  baselines,  a  linear  system  of  n 
equations  and,  principally  2  unknowns  (the  velocity  Vr  and  the  DEM-quantities)  is  established  and  a  joint 
estimation  of  both  DEM  errors  and  target  velocity  can  be  carried  out  on  a  pixel-by-pixel  basis.  In  order  to 
reduce  the  effect  of  atmospheric  patterns  and  compensate  for  baseline  errors,  the  low  frequency  components 
should  be  previously  removed  from  Oj.  The  processing  is  then  performed  on  a  pixel-by-pixel  basis  and  the 
final  result  is  a  velocity  field  of  the  area  of  interest  together  with  an  improved  DEM.  In  order  to  compensate 
for  possible  baseline  errors  and/or  very  low  frequency  spatial  phase  distortion  due  to  atmospheric  inhomoge¬ 
neities  the  linear  component  of  the  residual  phase  values  was  removed  from  all  the  available  data. 


Fig.  16  shows  the  principal  structure  of  an  overall  PSI  system.  (Adam,  et.  al.  2003). 
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Fig  16.  The  structure  of  the  overall  PS  system  of  DLR. 
Each  box  illustrates  a  sub-processing  system. 


The  meaning  and  the  tasks  of  the  single  Boxes  are  as  follows: 

The  Data  Import  box  starts  the  initial  data  setup  and  selects  the  area  of  interest.  The  InSAR  processing  box 
represents  the  normal  Interferometry  processing  procedure.  The  master  scene  is  selected  in  order  to  start  the 
InSAR  Processing.  Parameters  as  the  effective  baseline,  the  acquisition  date,  the  Doppler  Centroid 
frequency  and  the  season  of  the  acquisition  form  the  selection  criteria.  During  the  InSAR  processing  the 
observation  of  geometry  parameters  as  the  height-to-phase  conversion  factor,  the  flat-earth  phase,  the  range 
distances  and  the  look  angle  are  computed.  Furthermore,  the  scenes  are  co-registered  and  re-sampled  with 
respect  to  the  master  scene  and  the  Interferogrammes  are  generated.  Because  of  the  utilisation  of  point 
scatterers,  a  spectral  shift  filtering  is  not  necessary. 

In  the  Differential  Interferometric  Processing  module  the  observation  geometry  of  the  radar  acquisition  is 
simulated.  A  DEM  and  precise  orbits  are  used  as  input  in  order  to  correct  for  the  Interferogrammes.  Any 
available  DEM  can  be  used  to  reduce  the  topographic  signal.  Exemplary  the  DEM’s  generated  from  the 
X-SAR  shuttle  mission  can  be  used  ore  a  DEM  derived  from  the  ERS  tandem  mission.  This  module  also 
determines  the  co-registration  and,  therefore,  it  can  be  used  for  that  processing  step  also.  With  the  PS 
Interferometry,  the  co-registration  becomes  more  difficult  than  in  normal  Interferometry  considered  in 
Chap.4.  The  reason  is  the  long  time  separation  (up  to  ten  years)  between  the  acquisitions  and  the  immense 
temporal  de-correlation  in  non-urban  areas. 


The  do  calibration  of  the  scenes  has  to  be  introduced  for  the  PS  detection.  To  analyse  the  temporal  back- 
scattering  behaviour  of  the  point  scatterers  it  is  necessary  to  compare  their  intensity  by  the  calculation  of  the 
radar  cross  section.  This  module  corrects  for  the  processor  gain  constant,  the  antenna  pattern  and  the  range 
spreading  loss.  The  Permanent  Scatterer  Detection  and  identification  is  performed  in  the  co-registered 
calibrated  scenes  by  a  homogeneity  test  (Ferretti  et  al.,  2001).  Aim  is  to  find  as  many  scatterers  as  possible 
because  a  subsidence  pattern  and  the  atmosphere  have  to  be  sampled  spatially  as  dense  as  possible.  Flowever, 
unreliable  points  causing  incorrect  estimation  have  to  be  avoided.  Therefore,  the  temporal  analysis  of  the 
differential  phases  is  restricted  to  point  scatterers  with  a  high  SNR  and  a  long-time  stable  backscattering 
behaviour.  These  are  usually  man-made  features  (Usai  et  al.,  1999)  and  have  to  be  detected  in  the  scene.  The 
extraction  of  the  PS  information  reduces  the  amount  of  data  from  more  than  100  GB  to  less  than  300  MB. 
Range  and  azimuth  position  of  the  scatterers  in  the  scene,  the  calibrated  intensity  ,  the  look  angle,  the  Dop¬ 
pler  Centroid  frequency,  the  height-to-phase  conversion  factor,  the  temporal  baseline,  the  scatterer  sub  pixel 
position  and  the  differential  phase  in  a  compact  form  are  the  result  from  this  processing  step  which  can  be 
considered  as  a  raster  to  vector  converter. 
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Both  modules  Estimation  and  Visualisation:  The  measured  differential  phase  is  composed  of  contributions 
from  the  uncompensated  topography  Because  the  phase  is  measured  modulo  2n  a  non-linear  inversion  prob¬ 
lem  is  given.  It  is  solved  by  utilising  the  different  behaviour  of  the  contributions  regarding  the  acquisition 
parameters  effective  baseline,  temporal  baseline,  range  and  azimuth  location  of  the  scatterers.  The  visualisa¬ 
tion  box  visualises  the  different  spatial  properties  for  a  subsidence,  for  an  orbit  error,  for  the  atmosphere  and 
for  noise.  A  relative  estimation  between  point  scatterers  located  closely  to  each  other  reduces  the  influence  of 
the  atmosphere  and  of  orbit  errors.  Utilising  a  periodogram  a  constant  relative  subsidence  rate  can  be  esti¬ 
mated.  These  relative  estimates  are  transformed  into  a  global  subsidence  map  by  a  2-D  integration  procedure 
based  on  a  LS-adjustment. 

Fig.17  shows  a  section  of  a  radar  image  of  Berlin.  The  PS’s  are  marked  by  coloured  dots,  the  colours  repre¬ 
sent  the  subsidence  rate.  A  second  example  is  shown  in  Fig.  18.  The  area  around  the  Olympic  stadium  in 
Berlin  is  subject  to  ground  water  regulation  activities.  This  area  shows  up  as  red  PS’s,  i.e.,  the  ground  level  is 
lifting  by  about  3  mm/year.  During  the  lecture  on  applications  additional  examples  will  be  discussed. 


Permanent  Scatterer  Interferometry 
Subsidence  Rate 


±2ct  Ambiguities 


Fig.  17  Subsidence  rate  derived  by  the  PS  method.  The  marked  point  is  a  subway  station  and  settles  at  a  ve¬ 
locity  of  2.5  mm/year  compared  to  its  neighbourhood.  63  ERS  data  sets  from  1991  to  2000  have  been  used. 
The  plot  shows  the  relatively  low  noise  of  the  measurement.  The  ±  25  ambiguities  areplotted  for  easier  inter¬ 
pretation. 


Fig.  18  Inflation  due  to  ground  water  regulation  at  the  Berlin  Olympic  stadium.  PS  measurements  have  been 
spatially  interpolated  to  allow  for  a  better  visual  interpretation;  green:  no  subsidence  /  inflation,  red:  inflation 
of  4  mm/year,  cross:  reference  point,  assumed  to  be  stable:  (Kampes  et  al.,  2003). 

8.  Multifrequency  SAR  Interferometry 

A  multi-frequency  SAR  system  offers  a  great  potential  to  improve  the  unwrapping  performance  of  Interfero- 
grammes. 
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The  SIR-C/X-SAR  Shuttle  Mission  flown  twice  in  April  and  October  1994  provided  for  the  first  time  two 
pass  interferometric  data  in  X-,  C-,  and  L-band  gained  simultaneously  over  some  test  areas.  Respective  in¬ 
vestigations  have  shown  that  the  availability  of  multi  frequency  information  permits  a  drastic  improvement 
of  the  phase  unwrapping  performance  investigation  [Moreira  J.  et.  al  1995].  In  the  Fig  19  to  21  the  sets  of 
tree  co-registered  slant  range  images,  coherence  images,  and  a  set  of  slant  range  interferometric  phases  over¬ 
lapped  with  the  interferometric  magnitude  images  in  L-,  C-,  and  X-Band  is  depicted. 


Fig.  19  Set  of  co-registered  slant  range  SAR  images  in  L-,  C-,  and  X-band.  The  image  size  is  66.8  km  in  azi¬ 
muth  and  9  km  in  ground  range  (16000  radar  pulses  and  1024  range  bins  in  X-band  and  2048  range  bins  in 
L-  and  C-band). 


Fig.20  Set  of  co-registered  interferometric  Phase  fringes  in  L-,  C-,  and  X-Band  overlapped  with  the  respec¬ 
tive  magnitude  images  shown  in  Fig.  18. 
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Fig.21  Set  of  co-registered  slant  range  coherence  images  in  L-,  C-,  and  X-band.  White  indicates  highcoher- 
ence  and  black  low  coherence.  It  is  evident  that  1-Band  has  a  high  coherence  with  the  exception  of  certain 
small  areas 

Considering  an  image  pair  which  meets  the  baseline  and  Doppler  centroid  requirements,  the  interferometric 
phase  quality  is  limited  by  the  temporal  de-correlation  as  has  been  mentioned  already.  The  height  accuracy 
depends  on  both  the  wavelength  and  the  coherence.  For  areas  with  the  same  coherence  for  L-,  C-,  and  X- 
band,  the  height  derived  from  the  X-band  phase  is  approximately  two  times  more  accurate  than  C-band  and 
eight  times  more  than  L-band  due  to  the  wavelength  ratios. 

In  the  example  of  SIR-C/X-SAR,  the  basic  shape  of  the  fused  DEM,  i.e.  the  low  frequency  part,  is  mainly 
derived  from  L-band  due  to  its  high  coherence.  The  high  frequency  part  of  the  DEM  is  derived  from  C-  and 
X-band  due  to  their  small  wavelengths.  This  fusion  is  done  by  using  a  simple  Kalman  filter  and  takes  into 
account  the  respective  coherences. 

The  phase  gradient  estimation  can  be  improved  drastically  in  the  presence  of  critical  areas,  i.e.  low- 
coherence  areas,  if  an  a  priori  knowledge  of  the  terrain  topography  is  available.  The  main  idea  is  to  improve 
the  unwrapping  procedure  in  C-  and  X-band  by  using,  as  a  priori  information,  the  L-band  unwrapped  values. 
The  selection  of  the  L-band  data  is  related  to  its  higher  coherence  and  the  lower  sensitivity,  compared  to  C- 
and  X-band  data,  of  the  interferometric  phase  to  the  variations  of  the  height  profile. 

A  block  diagram  of  the  multi  frequency  unwrapping  procedure  is  depicted  in  Fig  22.  The  first  step  of  the 
proposed  procedure  is  represented  by  the  filtering  operation  of  the  L-band  wrapped  phase  (via  a  3  x  3  pixel 
pivoting  median  filter)  followed  by  the  phase  unwrapping  step.  The  L-band  unwrapped  data  are  used  to  im¬ 
prove  the  C-and  X-band  unwrapping  results  following  the  procedure  depicted  in  Fig.22. 
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Er.  Xx/Ax 


Fig.22  Block  diagram  of  the  multi-frequency  phase  unwrapping  procedure.  (<cpL  >,  (<(pC  >  and  <  cpX  >  are 
the  wrapped  phases  of  the  L-,  C-,  and  X-band,  respectively,  without  the  flat-earth  component;  ,A.L,  Xc,  Xx  are 
the  wavelengths.  Ef  is  the  flat-earth  function  scaled  for  X-band  also  used  in  the  L-  and  C-band  with  the 
respective  wavelength  scaling.  *  denotes  the  phase  difference  calculated  in  the  complex  domain.  ®L,  ®c,®x 
are  the  unwrapped  original  phases  in  the  respective  Bands. 

First  of  all,  the  resulting  L-band  phase  is  scaled  with  respect  to  the  C-band  wavelength  and  subtracted  from 
the  C-band  measured  phase.  This  operation  is  carried  out  in  the  complex  domain  so  that: 

•  both  rescaled  L-band  and  C-band  phases  are  transformed  into  complex  signals  by  using  an  unitary 
amplitude, 

•  the  difference  operation  is  realized  by  multiplying  the  complex  C-band  signal  by  the  complex 
conjugate  of  the  L-band  signal,  and 

•  the  wrapped  phase  difference  is  obtained  from  the  arctangent  of  the  above-mentioned  complex 
multiplication 

The  result  of  the  difference  operation  is  filtered  with  a  3  x  3  pivoting  median  filter  and  then  unwrapped.  It 
should  be  noted  explicitly  that  the  pivoting  median  filter  will  benefit  from  the  improved  phase  field 
stationary  due  to  the  performed  subtraction.  Finally,  the  computed  phase  is  added  to  the  L-band  reference 
signal. 

The  same  strategy  is  used  to  unwrap  the  X-band  phase  Interferogramme.  The  only  difference  in  this  case  is 
that  the  previously  computed  C-band  phase,  instead  of  the  L-band  one,  is  used  as  a  reference  pattern.  The 
overall  procedure  block  diagram  is  depicted  in  Fig.  23. 
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Fig.  23  Block-diagram  of  the  Phase  Fusion  Algorithm 
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Fig.  24  Unwrapped  phase  profiles(  azimuth  direction)  filtering  without  unwrapped  phase  profiles  of  L-,  C-, 
and  X-band  in  azimuth  direction  applying  multi-frequency  phase  unwrapping  without  filter  before  phase 
unwrapping.  The  L-  and  C-band  phases  were  scaled  with  respect  to  the  X-band  wavelength  for  comparison 
purposes.  The  x-axis  gives  the  azimuth  dimension  in  pixels.  Each  pixel  has  a  size  of  20.9  m  (3200  pixels 
correspond  66.864  km).  The  y-axis  unit  is  radians. 
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Fig.  25  Unwrapped  phase  of  L-,  and  X-band  profiles  in  azimuth  direction  applying  multi-frequency  phase 
unwrapping  with  the  cascade  of  two  filters  with  a  window  of  5x5  before  phase  unwrapping.  The  L-  phase 
was  scaled  with  respect  to  the  X-band  wavelength  for  comparison  purposes.  The  x-axis  gives  the  azimuth 
dimension  in  pixels.  Each  pixel  has  a  size  of  20.9  m  (3200  pixels  correspond  66.864  km).  The  y-axis  unit  is 
radians. 


Interferometric  SAR  DEM 


Test  site:  Mt.  Etna  /  Sicily 


Fig  26  Perspective  view  of  the  geo-coded  DEM  of  the  Mount  Etna  area.  The  size  of  the  x-y  plane  is  25  km  x 
25  km,  the  unit  along  the  z-axis  is  meters. 

To  demonstrate  the  improvement  achieved  by  using  the  multi- frequency  information  the  corresponding  re¬ 
sults  obtained  by  independently  unwrapping  the  three  Interferogrammes  without  any  filtering  and  filtered  are 
depicted  in  Fig.24  and  Fig.  25  which  show  slices  through  an  unwrapped  Interferogamme  before  and  after 
filtering  (next  page).  The  comparison  makes  evident  the  joint  effects  of  the  filtering  step  combined  with  the 
introduction  of  the  multi- frequency  information.  In  Fig.  24  there  are  still  some  reconstruction  problems 
caused  by  extremely  low-coherence  areas  present  in  the  X-band  data  where  the  bias  could  not  be  completely 
avoided.  A  better  result  can  be  obtained  by  applying  stronger  data  filtering.  This  is  shown  in  Fig.  25,  where 
the  L-band  reference  is  compared  to  a  highly  filtered  X-band  unwrapped  pattern  (retrieved  in  the  multi- 
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frequency  phase  unwrapping  procedure  by  using  the  cascade  of  two  5x5  median  filters).  Obviously,  the  im¬ 
provement  obtained  is  at  the  expense  of  a  loss  in  the  height  resolution  in  X-band.  In  Fig.  26  is  the  resulting 
DEM  of  the  Mount  Etna  depicted. 


9.  SAR  Tomography 

Interferometry  determines  by  means  of  the  phase  differences  between  two  different  sensor  positions  the 
terrain  topography.  Although  it  is  possible  in  this  way  to  find  a  three-dimensional  (3-D)  surface 
representation,  the  distribution  of  the  different  scatterers  in  the  height  direction  at  a  fixed  range  and  azimuth 
position  remains  unknown.  Contrary  to  this,  tomographic  techniques  enable  a  real  geometric  resolution 
capability  in  the  height  direction  and  introduce  new  possibilities  for  many  applications  and  inversion 
problems.  Even  misinterpretations  in  SAR  images  caused  by  layover  and  foreshortening  effects  can  be 
solved  by  the  tomographic  processing. 

The  combination  of  interferometric  and  polarimetric  techniques  enables  the  separation  of  different  scattering 
mechanisms  within  a  resolution  cell  and  at  the  same  time,  the  estimation  of  the  associated  heights.  The  phase 
decomposition  performed  in  polarimetric  SAR  interferometry  uses  the  full  polarimetric  data  and  estimates 
for  each  image  pixel  three  orthogonal  scattering  mechanisms.  For  example,  for  a  vegetated  area  the  three 
different  scattering  mechanisms  can  be  described  as  (in  a  simplified  model): 

1)  a  single  bounce  contribution  from  the  ground; 

2)  a  double  bounce  contribution  from  the  trunks;  and 

3)  a  diffuse  (volume)  contribution  from  the  tree  leaves  and  branches. 

However,  it  is  not  possible  with  this  technique  to  separate  the  contributions  of  the  same  scattering  mecha¬ 
nism  distributed  over  different  heights. 

The  basic  idea  of  Tomography  is  to  form  a  second  synthetic  aperture  in  the  direction  of  the  translational 
movement  (see  Fig.  27).  This  direction,  denoted  as  normal  direction  n,  is  perpendicular  to  the  line-of-sight 
(range  or  across-track)  R  and  azimuth  (flight)  directions  x. 
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Fig.  27  Tomography  for  Scattering  Analysis  (Reigber,  A.,  1999)  Left  Part:  Top  Principle  of  SAR  Tomogra¬ 
phy  with  Multiple  horizontal  Flights  as  Interferometry  with  Several  Vertical  Baselines,  below  real  flight 
passes  of  a  Do-228.  Right  part:  E-SAR  results.  Top:  3-dim  forest.  Middle  SAR  image,  yellow  line  indicates 
the  flight  path.  Below:  Colour  coded  different  scattering  mechanisms  are  obtained  with  polarization.The 
buildings  roof  is  bright,  the  blue  single  bounce  is  the  comer  reflector,  red  double  bounces  happen  at  parts  of 
the  buildings  roof  and  at  the  lower  parts  of  trunks,  in  tree  crowns  and  under  story  volume  green  scattering 
diffuse  crosspol  scattering  dominates.  Very  below  is  artistic  view  of  the  imaged  area. 


The  synthetic  aperture  in  azimuth  direction  is  associated  with  a  variation  of  the  so-called  squint  angle,  while 
the  synthetic  aperture  in  the  normal  direction  is  associated  with  a  variation  of  the  off-nadir  angle  for  each 
object..  Therefore,  Tomograhy  can  be  interpreted  also  as  a  multi  baseline  Interferometry  with  the  number  of 
tracks  being  much  greater  than  two. 


Remarkable  is  the  deviation  of  the  real  flight  passes  in  Fig.  27  from  ideal  strait  lines.  In  connection  with  the 
former  considerations  on  the  influence  of  geometrical  errors  in  Interferometry  arrangements  depicted  in  Ta¬ 
ble  3  for  example  this  points  out  that  motion  compensation  procedures  are  indispensable  for  exact  inter- 
feometric  measurements.  This  will  be  discussed  in  the  Systems  Chapter  again. 

10.  SAR  Transponder  Techniques  for  Absolute  Altitude  Calibration 

Calibration  for  interferometric  SAR  systems  is  a  conditio  sine  qua  non  as  it  is  for  all  other  SAR  SAR  sytems. 
Mainly,  in  that  context  the  calibration  of  the  amplitude  of  SAR  signals  is  considered.  An  absolute  altitude 
calibration  in  interferometric  SAR  can  be  performed  via  Transponder  Techniques.  The  use  of  both  active  and 
passive  transponder  for  SAR  Calibration  is  a  common  and  well  developed  techniqueThis  will  be  discussed  in 
more  detail  in  the  Systems  chapter.  But  the  absolute  height  calibration  is  more  difficult  than  it  seems  under  a 
first  look. 

In  spite  of  the  high  resolution  achieved  by  SAR,  the  resolving  power  is  often  insufficient  to  detect,  identify 
and  classify  small  objects  in  hard  clutter  environment.  In  such  environments  transponders  with  special  cod¬ 
ing  can  be  used  to  identify  objects  in  SAR  images  via  tagging.  The  transponder  modulates  the  signal  re¬ 
ceived  from  the  satellite  with  a  characteristic  code  in  range  as  well  as  in  azimuth  [Hounam.  D,  1996]  using 
signal  modulation  techniques,  which  affect  the  SAR  signal  in  both  range  and  azimuth  and  decouple  the 
transponder  signal  from  the  background  signal.  This  causes  the  transponder  signal  to  be  defocused  and 
strongly  suppressed  in  the  SAR  image.  With  a  separate  processing  the  transponder  and  the  object  on  which  it 
is  mounted  can  thus  be  detected,  located  and  identified  even  when  situated  in  a  high  clutter  area.  Each  trans- 
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ponder  can  be  given  an  individual  code.  In  principle,  the  transponder  can  be  very  small  and  lightweight.  This 
offers  advantages  for  both  radiometric  and  three-dimensional  geometric  calibration.  The  retransmitted  signal 
can  be  coded  with  respect  to  the  transponder’s  three-dimensional  exact  position.  This  allows  an  exact  deter¬ 
mination  of  its  three  dimensional  position  especially  of  its  absolute  altitude  also.  (If  the  transponder  is 
equipped  additionally  with  a  GPS  receiver  chip  the  GPS  information  can  be  inserted  in  the  signal  transmitted 
to  the  SAR  also  which  would  allow  an  unattended  position  determination).  [Hounam,  D.  et.  al.,  1998,  Bloet- 
scher,  H.  et.al.  1998]. 

In  the  upper  part  of  Fig  28  a  Transponder  signal  within  a  highly  cluttered  ERS-1  scene  is  depicted  in  both 
situations  before  (left)  and  after  (right)  decoupling  of  the  signal  from  the  imaged  background.  The  strong 
suppression  of  the  transponder  in  the  SAR  Image  is  evident  as  well  as  the  strong  evidence  of  the  Trans¬ 
ponder  signal  after  Processing  which  takes  into  account  the  special  Transponder  code.  The  upper  images  are 
the  three  dimensional  SAR  signals  of  the  ERS-1,  the  lower  image  is  a  conventional  SAR  imge  of  the  same 
scene.  The  small  white  circle  in  the  middle  of  the  lower  images  gives  the  Position  of  the  Transponder. 


Fig.  28  Really  measured  example  of  a  coded  Transponder  in  a  clutter  surrounding  of  an  ERS-1  scene.  Upper 
left  and  lower  image:  the  Trans-ponder  is  shadowed  completely  and  not  detectable.  Upper  right  image:  The 
Transponder  clearly  shows  up  after  a  specially  converted  processing. 
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